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Abstract 

It has been shown recently that phantom-like effect can be realized on the normal 
branch of the DGP setup without introduction of any phantom matter neither in the 
bulk nor on the brane and therefore without violation of the null energy condition. It 
has been shown also that inclusion of the Gauss-Bonnet term in the bulk action modifies 
this picture via curvature effects. Here, based on the Lue-Starkman conjecture on the 
dynamical screening of the brane cosmological constant in the DGP setup, we extend 
this proposal to a general DGP-inspired f(R, <p) model that stringy effects in the ultra- 
violet sector of the theory are taken into account by inclusion of the Gauss-Bonnet 
term in the bulk action. We study cosmological dynamics of this setup, especially its 
phantom-like behavior and possible crossing of the phantom divide line especially with 
a non-minimally coupled quintessence field on the brane. In this setup, scalar field and 
curvature quintessence are treated in a unified framework. 
PACS: 98.80.-k, 95.36,+x, 98.80.Cq 
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1 introduction 



Recent evidences from supernova searches data [1,2], cosmic microwave background (CMB) 
results [3-5] and also Wilkinson Microwave Anisotropy Probe (WMAP) data [6,7,8], show 
an positively accelerating phase of the cosmic expansion today and this feature shows that 
the simple picture of the universe consisting of the pressureless fluid is not enough to de- 
scribe the cosmological dynamics. In this regard, the universe may contain some sort of the 
additional negative-pressure dark energy. Analysis of the WMAP data [9-11] shows that 
there is no indication for any significant deviations from Gaussianity and adiabaticity of the 
CMB power spectrum and therefore suggests that the universe is spatially flat to within the 
limits of observational accuracy. Further, the combined analysis of the WMAP data with 
the supernova Legacy survey (SNLS) [9], constrains the equation of state u^ e , corresponding 
to almost 74% contribution of dark energy in the currently accelerating universe, to be very 
close to that of the cosmological constant value. In this respect, a ACDM ( Cosmological 
constant plus Cold Dark Matter) model has maximum agreement with the recent data, but 
with an unnaturally small and fine-tuned value of A [4]. From the general relativity view- 
point, such models simply modify the right hand side of the Einstein's field equations by 
including in the stress-energy tensor something more than the usual matter and radiation 
components. As a radically different approach, one can also try to leave the source side un- 
changed, but rather modifying the geometric sector in the left-hand side. In this viewpoint, 
the cosmic speed up can be interpreted as a first signal of the breakdown of the laws of 
physics as encoded by the standard General Relativity. This is essentially reasonable: the 
General Relativity has been experimentally tested only up to the Solar System scale and 
there is no a priori theoretical motivation to extend its validity to extraordinarily larger 
scales such as the cosmological ones. Extending General Relativity without giving up its 
positive results, opens the way to a large class of alternative theories of gravity ranging from 
extra-dimensions [12] to non-minimally coupled scalar fields [13,14,15]. In particular, fourth 
order theories [16,17,18] based on replacing the scalar curvature R in the Einstein-Hilbert 
action with a generic analytic function f(R) have attracted a lot of attention in the last 
decade. Also referred to as f(R) gravity, these models have been shown to be able to both 
fit the cosmological data and evade the Solar System constraints in several physically inter- 
esting cases [19-28]. In this sense, modified gravity provides the very natural gravitational 
alternative for dark energy. Also modified gravity presents very natural unification of the 
early-time inflation and late-time acceleration due to different role of gravitational terms rel- 
evant at small and at large curvature [29,30,31]. Moreover, some models of modified gravity 
are predicted by string/M-theory considerations [16,17,18]. 

In the framework of braneworld scenarios, a possible approach to realize the late-time 
acceleration of the universe is to consider modification of gravity on the appropriate scales. 
This can be achieved for instance by a weakening of gravitational interaction at large scales. 
Weakening of gravity on large scales provides an effective negative pressure that drives the 
late-time acceleration of the universe. The Dvali, Gabadadze and Porrati (DGP) scenario 
[32] is one of the simplest scenario in this regard. In this model, our universe is a brane 
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embedded in a Minkowski 5D bulk space-time. Standard matter and gauge interactions are 
trapped on the brane and only gravity experiences the full bulk. Cosmological dynamics in 
this setup follows two distinct branches of solutions: self-accelerating and the normal branch 
[33,34,35]. Late-time acceleration can be explained in the self-accelerating branch without 
invoking any unknown dark energy component. However, the normal branch requires a 
dark energy component to accommodate the current observations [36-45]. Although the 
self-accelerating branch explains the late-time acceleration in a fascinating manner, it suffers 
from serious theoretical problems such as the ghost issue [46]. Fortunately the normal branch 
is free from instabilities. Recently it has been shown that the normal branch of the DGP 
braneworld scenario which is not self-accelerating can realize phantom-like behavior without 
introducing any phantom field on the brane [36-45]. By phantom-like behavior we mean 
an effective energy density which is positive and grows with time and its equation of state 
parameter stays always less than —1. Other extensions to realize this phantom-like behavior 
are studied in Refs. [47-66]. On the other hand, dark energy models with non-minimally 
coupled scalar field and other extensions of scalar-tensor theories have been studied widely 
some of which can be found in references [67-71]. 

In the spirit of scalar-tensor dark energy models, our motivation here is to show that a 
general f(R, 4>) GBIG-inspired scenario can account for late-time acceleration and crossing 
of the phantom divide line in some suitable domains of model parameters space. To show 
this feature, first we study cosmological dynamics of f(R, <j>) GBIG-inspired scenario briefly. 
In the minimal case, motivated by modified theories of gravity, some authors have included 
a term of the type l R n in the action [72-77]. This extension for some values of n has the 
capability to explain late-time acceleration of the universe in a simple manner. It is then 
natural to extend this scenario to more general embeddings of DGP inspired scenarios. 

The purpose of this paper is to perform such a generalization to study both late-time 
acceleration and possible crossing of the phantom divide line in this setup, we construct 
a braneworld scenario with modified induced gravity on the brane and the Gauss-Bonnet 
term in the bulk action. We study cosmological dynamics of this setup and we explore 
phantom-like behavior in the normal branch of the model. We focus mainly on the potential 
of the type V(<f>) = \(p n for non-minimally coupled scalar field. Our setup for minimal 
case predicts a power-law acceleration supporting observed late-time acceleration. In the 
non-minimal case, by a suitable choice of non-minimal coupling and scalar field potential, 
one obtains accelerated expansion in some specific regions of parameters space. While a 
single minimally coupled scalar field in four dimensions cannot reproduce a crossing of the 
phantom divide line for any scalar field potential [47-66], a non-minimally coupled scalar 
field account for such a crossing [78]. In GBIG model, equation of state parameter of dark 
energy never crosses cu(z) = —1 line, and universe eventually turns out to be de Sitter phase. 
Nevertheless, in this setup if we include a single scalar field (ordinary or phantom) on the 
brane, we can show that equation of state parameter of dark energy can cross the phantom 
divide line [79] . Crossing of the phantom divide line with non-minimally coupled scalar field 
on the warped DGP brane has been studied recently [80]. Here we present an extension of 
these non-minimal dark energy models within GBIG-inspired f(R, (f>) braneworld scenario. 
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We use a prime for differentiation with respect to R. An overdot marks differentiation 
with respect to the brane time coordinate. We show that this model provides a natural 
candidate for explanation of several interesting aspects of the late-time acceleration without 
suffering from issues such as the ghost instability. Among these features, natural realization 
of the phantom divide line crossing is worth to noting. 



2 The Modified GBIG Scenario 

The GBIG gravity is a generalized braneworld scenario which contains both UV (ultra-violet) 
and IR (infra-red) limits in a unified manner: It contains stringy effect via Gauss-Bonnet 
(GB) term in the bulk action as the UV sector of the theory and Induced Gravity (IG) 
effect which becomes important in the IR limit [81-89]. Here we add a new ingredient to 
this scenario: we consider possible modification of the induced gravity on the brane in the 
general framework of f(R, (f>) gravity. The action of our modified GBIG scenario is written 
as follows 



S = -L [ d 5 xyJ-g( 5 )\R® -2A 5 + a 



R^ 2 - AR^R^ ab + R a 5 b \ d R {5)abcd 



+ / d A Xy^g 

J brane 



1 ^ 



1 



K 



2 K + ^f(R,cP)-X + L r , 



v 5 z,n, 4 

where a is the GB coupling constant. We assume the brane is located at y — 0, where y is 
the coordinate of the fifth dimension. K is the trace of the mean extrinsic curvature of the 
brane defined as follows 
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(2) 



and its presence guarantees the correct matching conditions across the brane. We denote 
the matter field Lagrangian by L m with the following energy-momentum tensor 



T, 



2 J^ + 9lluLr 



(3) 



The field equations resulting from the action (1) are given as follows 

Ml{Gi + a Hi)=hs8i + 8ir B T». 



(4) 



We note that a contribution of the form 



,(5) 



is inherent in the definition of t$ on the right. 



The corrections to the Einstein Field Equations originating in the GB term are represented 
by the Lovelock tensor [90] 
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The energy-momentum tensor localized on the brane is 



+M 2 [D»D v f'(R, 0) - 5»D a DJ'(R, </>)] 



A 



(6) 

where -D M is the covariant derivative with respect to . The corresponding junction 
conditions relating quantities on the brane follows [91,92] 
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(7) 



To formulate cosmological dynamics on the brane, we assume the following line element 

ds 2 = -n 2 (y, t)dt 2 + a 2 (y, t)j ij dx i dx j + b 2 (y, t)dy 2 (8) 

where 7^ is a maximally symmetric 3-dimensional metric defined as 7^ = 5^ + k ^'^2 where 
k = —1, 0, +1 parameterizes the spatial curvature and r 2 = x^x 1 . y is the coordinate of extra 
dimension and the brane is located at y — 0. The junction conditions on the brane now 
gives the following expressions 
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If we choose a Gaussian normal coordinate system so that b 2 (y, t) = 1, these equations with 
non- vanishing components of the Einsteins tensor in the bulk yield the following generaliza- 
tion of the Friedmann equation for cosmological dynamics on the brane ( see for instance 
[93-98] for the general machinery of this derivation ) 
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where p tot = p + p c + p</>- The bulk contains a black hole mass and a cosmological constant 
so that $ is defined as follows 
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where T is an integration constant. We note that r c is the crossover scale in the DGP model 
which measures the strength of the induced gravity effect on the brane and is related to the 
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four and five dimensional gravitational constants by r c = Also, the parameter a which 
is positive measures the strength of the GB curvature effect on the brane and has dimension 
of length square. If a = 0, the model reduces to the DGP model with modified induced 
gravity, while for r c = we recover the Gauss-Bonnet braneworld model. Here we restrict 
our study to the case where the bulk black hole mass vanishes (that is, T = 0) and therefore 
$ + 2a$ 2 = In this case the bulk cosmo logical constant is given by A5 = ^ + ^pr, where 
I is the bulk curvature. Assuming A5 = corresponding to a Minkowski bulk, for a spatially 
flat FRW brane (K = 0), the Friedmann equation is given by 
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where A is the brane tension. We restrict ourselves to the normal branch of this DGP-inspired 
model. The energy density on the brane is given by a CDM component ( with energy density 
p m ) and a cosmological constant A with p = p m +A and a term p c originating on the modified 
induced gravity as 
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where p m = p mo (l + z) 3 . With A 5 = we find $ = and $ = Our forthcoming 

arguments are based on the choice $ = 0. Now, if we define the cosmological parameters as 

2 

and Qtf, = then the 
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Friedmann equation on the brane can be expressed as follows 
1 



E\z) 



f'(RA) 



2Jn rc E(z)[i + n a E 2 ( z )} + n m (i + zf + + n c + ^ A 



(18) 



where E(z) = j^. Evaluating the Friedmann equation at z = gives a constraint equation 
on the cosmological parameters of the model as follows 



1 + 2Jn rc (i + n a ) = tt m + n c + n A + 



(19) 
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which can be written formally as 



i + 2^/n rc (i + n a ) = n m + n tot . (20) 

We note that a contribution to this equation from f'(R, (f>)\ z =o nas been normalized to unity. 
This is reasonable since we finally have to choose an ansatz for f(R, 4>) where R and <p are 
related to the redshift via scale factor and its time derivative. This quantity computed at 
z = gives a number which we normalize it to unity. This constraint equation implies that 
Q m + Q c + Q\ + < 1 is unphysical in the model parameter space. By imposing the 
condition that the universe is currently accelerating, the deceleration parameter defined as 
q = —[jj2 + 1] takes the following form 



f'(R sh\ J- 3^m-Hp(l+z) 3 _ «t(Pc+P0) 

J e(z) m-EHz) 



H [2f'(R, ^)E{z) + 2^(1 + Q a E^(z)) + A^Kn a E\z)} 

To proceed further, we should specify the functional form of f(R,(f>). One of the most 
pressing problems of f(R, 0) theories is the need to escape the severe constraints imposed 
by the Solar System tests [23], [99-101]. The Hu-Sawicki model with [25] 

«o(i)" 



f(R)-R-R Ci + Po& 

is among those models that have good agreement with Solar-System tests. In which follows, 
we choose the ansatz 

/(^)=(fl- 0'|I my (22) 

where the independent part is from Hu-Sawicki scenario and in the ^-dependent part we 
have included a non-minimal coupling term. We adopt also the ansatz <p(t) = <pot~^ , where 
jj, is a positive constant and we set n = 0.9. For scale factor we adopt the following ansatz 

a(t) = (e + ^-)^ (23) 

where that v ^ 1 [102]. Figure 1 shows the behavior of q versus the redshift in this model. 
The universe enters the accelerating phase at z ~ 0.44. 

The modified Raychaudhuri equation follows easily from the Friedmann equation of the 
brane and the conservation of the brane energy density 

H f'(R, 4>)E 2 {z) + 3Q m ff (l + zfE{z) - 

— 6 Q (od) 

Hi H [2f\R,<p)E(z) + 2^- c (l + n a E^z))+A^:n a E\z)Y K ' 

Figure 2 shows the variation of with redshift. It is always negative and therefore there is 

no super-acceleration and big-rip singularity in this braneworld scenario. Since H < 0, the 
Hubble parameter decreases as the brane expands. 
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3 The Phantom-like Behavior 



We will show in this section that a phantom-like behavior can be realized on the brane 
and this occurs without including any exotic matter that violates the null energy condition 
not on the brane nor in the bulk. The phantom-like behavior is based on the definition of 
an effective energy density which is corresponding to a balance between the cosmological 
constant and geometrical effects encoded on the Hubble rate evolution [36-45], [4 7-66]. The 
phantom-like behavior on the brane is based on the definition of an effective energy density 
p e ff which increases with cosmic time and an effective equation of state parameter always 
less than —1, u) e ff < —1. More precisely, the effective description is inspired in writing down 
the modified Friedmann equation of the brane as the usual relativistic Friedmann equation 
so that [47-66] 



+ Peff) 



(25) 



Using equations (13) and (25), we find 



Peff 



f'(R, 



2Jn rc E( z )(i + n a E\z)) + n m (i + z y + n c + ^ + n A - n m (i + z) 



.(26) 



By definition, phantom-like prescription breaks down if p e ff < 0. Figure 3 shows the vari- 
ation of p e ff versus z in our proposed setup. It is always positive and grows with cosmic 




Figure 3: p e ff versus redshift. 



time, therefore effective energy density shows the phantom-like behavior. In this respect, the 
important feature of this scenario is the fact that this phantom-like behavior never breaks 
down at least in the proposed subspace of the model parameter space. Nevertheless, there 
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are other subspaces of the model parameter space that phantom-like prescription may be 
broken in these subspaces in the same way as has been mentioned by Bouhamdi-Lopez in Ref. 
[47-66]. The main point here is the fact that it is possible essentially to save phantom-like 
prescription for a wide range of redshifts. 

The effective equation of state parameter u e ff, can be defined using the conservation 
equation of the effective energy density on the brane 

p eff + 3H(l + uj eff )p eff = 0. (27) 

The effective energy density evolve as follows 

[2HH ^~(1 + Q. a E\z)) + 2Q a E 2 (z)HH + 3Q m E(z)(l + zfH* - <£f£±] 



Peff = ~2 



+ ( /W ))2 ^ m (l + z) E(z)H 



.(28) 
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_ f(R, 0)// o 2 [-2 v ^:(i + n a E 2 (z)) + n m (i + zf + n c + ^ + n A ] 



3Q m (l + z) 3 E(z)Hl,(29) 



we find 



_ 3 j2ijf/ 0v /^(l+n a E 2 ( Z ))+2a a E 2 ( 2 )ijHo+3n m £(2)(l+ Z ) 3 g3- ^ Pe 3 +p -» ] 

f'(R,<i>) + y\ 



1 + ^ = wfrk • (30) 

Figure 4 shows the plot of 1 + u e ff versus the redshift. The universe enters to the phantom 
phase in z ~ 0.25. As this figure shows, the phantom-like behavior has no break down 
at small redshifts in contrast to the ADGP-GB model [47-66]. So, our model realizes a 
phantom-like behavior without introducing any phantom matter on the brane or in the bulk. 
It is important to note that this setup explains also a smooth transition to phantom phase, 
the so called phantom- divide line crossing. Our main achievement here is the realization 
of the phantom-like behavior without introducing any phantom matter not on the brane 
nor in the bulk. Nevertheless, the phantom-like behavior realized in this gravitational way 
should respect the null energy condition, at least in some subspaces of the model parameter 
space. So, we should check the status of the null energy condition in this setup. The 
null energy condition is respected if p e ff + p e ff > 0. Figure 5 shows a plot of p e ff + p e ff 
versus the redshift. As this figure shows, the phantom-like behavior in our proposed model 
occurs without violation of the null energy condition at least in some subspaces of the model 
parameter space. Especially, this condition is respected in the phantom-like region of the 
model. We emphasize that validity of the null energy condition in the phantom-like region of 
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Figure 4: 1 + w e ff versus the redshift. 

the scenario is important for our purposes in this paper. This is because with phantom fields 
usually the null energy condition breaks down. Here we realized a phantom-like behavior in 
a braneworld setup without violation of the null energy condition in the phantom-like region. 
It is important to note that in the non-phantom region, the situation is different and the 
validity of null energy condition is guaranteed by positivity of the effective energy density in 
this region. 



p+p 




z 



Figure 5: The status of the null energy condition in this setup: p e ft +p e ff versus the redshift. 



4 Lue-Starkman Screening of the Brane Cosmological 
Constant 

Lue and Starkman have shown that one can realize phantom-like effect, that is, increasing of 
the effective dark energy density with redshift, in the normal branch of the DGP cosmological 
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solution without introducing any phantom field [103] . The normal branch of the DGP model 
which cannot explain the self-acceleration, has the key property that brane is extrinsically 
curved so that shortcuts through the bulk allow gravity to screen the effects of the brane 
energy- momentum contents at Hubble parameters H ~ r^ 1 where r c is the crossover distance 
[103]. Since in this case H(z) is a decreasing function of the redshift, the effective dark energy 
component is increasing with redshift and therefore we observe a phantom-like behavior 
without introducing any phantom matter that violates the null energy condition and suffers 
from several theoretical problems. In the first step, in this section we study the phantom-like 
effect in the normal branch of a GBIG scenario. In other words, here we suppose that there 
is a quintessence field non-minimally coupled to modified induced gravity on the brane in 
addition to a Gauss-Bonnet term in the bulk action. In fact we incorporate the possibility of 
modification of the induced gravity on the brane in the spirit of f(R) gravity. We emphasize 
also that we have included a canonical (quintessence) scalar field to incorporate possible 
coupling of the gravity and scalar degrees of freedom on the brane. This provides a wider 
parameter space with capability to handle the problem more completely. In fact, inclusion 
of this quintessence field brings the theory to realize crossing of the phantom divide line 
[36-40]. Considering the normal branch of the equation (13) , we have 

H = -^-{pm + + pc) + (l + -aH), (31) 

where G = G e ff = (8nM%F'(R,(f))j and a prime denotes differentiation with respect to 
R. Now, we rewrite this equation in the following form 

H 2 = ^{p m + P4> + Pc) + ^(p^)- (32) 
By comparing this equation with the Friedmann equation (31), we find 

SnG ff A [f{R^\^H 8 2 

— (Pde) =3 ~ (1 + g-"# )• (33) 

As a result, in this GBIG-inspired model the crossover scale takes the following form 

f = 2Mf = * ^ ' ^ ' ^ ' 

Now as an enlightening example, we set for instance 

ry ( -EL) n 1 

m, 4>) = (r - \ + °:% )n ) (2 (i - ^ 2 )) (35) 



+ PoyR- c 

where the f(R) part is the Hu-Sawicki model [25]. With this choice, one recovers the general 
relativity for n — 0. For n 7^ 0, we obtain 

-^{p e iL) = I - '° (fle) r (1 + , (36) 
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where for a spatially flat FRW geometry, the Riici scalar is given by 



i? = 6- + 6(-) 2 . (37) 
a a 

Figure 6 shows the variation of p®^j versus n and the redshift in this GBlG-inspired F(R, (ft) 
model. As we see, the phantom-like behavior can be realized for all values of n constraint 
in the Hu-Sawicki model. In other words, for all values of acceptable n in the Hu-Sawicki 
model, the effective dark energy in this GBIG- inspired F(R, cf)) model has the phantom-like 
behavior. Figure 7 shows the variation of the effective dark energy density versus z and 
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Figure 6: Variation of the effective dark energy density versus n and z in a GBIG-inspired F(R, < 



model with F(R, <p) 
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the non-minimal coupling. By increasing the values of £, the effective dark energy density 
reduces, and for a fixed value of £ there is phantom-like effect for appropriate values of 
n. Also, figure 8 shows the variation of the effective dark energy density versus n and the 
non-minimal coupling. The phantom-like effect ( increasing the values of the effective dark 
energy density with cosmic time) can be realized for suitable ranges of n and £. 

As an important especial case, for F(R, (ft) = f(R) we find the screening effect in a 
general GBIG /(i?) -gravity 
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Figure 7: Variation of the effective dark energy density versus z and the non-minimal coupling. 



rho 




Figure 8: Variation of the effective dark energy density versus n and the non-minimal coupling. 
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Figure 9 shows the variation of the effective dark energy density versus n. There is phantom- 
like effect for all possible values of n in the Hu-Sawicki model. Figure 10 shows the variation 




Figure 9: variation of the pde versus n. 

of the effective dark energy density versus the redshift. The main result of these numerical 
solutions is that this modified GBIG scenario has the capability to realize the phantom-like 
prescription without introducing any phantom field neither on the brane nor in the bulk. 
We note that we have included a quintessence field non-minimally coupled to the modified 
induced gravity on the brane to have a wider parameter space and to include all possible 
situations. 

5 Dynamical Screening and Phantom-Like Mimicry with 
non-minimal coupling and curvature effect 

Our setup in this section contains a DGP-inspired braneworld scenario with curvature effect 
through the Gauss-Bonnet term in the bulk action and a non-minimally coupled scalar field 
on the brane. The ordinary matter localized on the brane is assumed to be a quintessence 
scalar field non-minimally coupled to the induced gravity on the brane. This canonical 
scalar field respects the null energy condition and we assume there is no phantom matter 
present on the brane or in the bulk. By virtue of dynamical screening, we show that this 
model accounts for a phantom-like behavior without violating the null energy condition. In 
addition, we show that in contrast to existing literature which argued that this phantom-like 
effect breaks down in some small values of redshift, there is no break-down of this behavior in 
the presence of the non-minimal coupling of the quintessence field . The action of our model 
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Figure 10: : pde versus the redshift. 



with Gauss-Bonnet and induced gravity effects and a scalar field non-minimally coupled to 
induced gravity on the brane is given as follows [104] (see also [81-84]) 



S = ^( d 5 x^gV) - 2A 5 + a(R^ 2 - 4R^R^ ab + R { H d R^ abcd 



+ 



2«§ 



J d A x^(f3(<j))R - 2Klg^V^VA - 4/^(0) 



(40) 



Jj/=0 



where a and r c are the GB coupling constant and IG cross-over scale respectively. We 
assume the brane is located at y = where y is the coordinate of the extra dimension. Also 
/3 '(</>) is a general non- minimal coupling of the scalar field and induced gravity on the brane. 
The cosmological dynamics of the model is given by the following generalized Friedmann 
equation [104] 



1 8 (tt2 $ K \ 



H 2 



c mH 2 + r c m- -M p+P4>+ x) 



.(41) 



This equation describes the cosmological evolution on the brane with tension and a non- 
minimally coupled scalar field localized on the brane. As usual, the bulk contains a black 
hole mass and a cosmological constant, so that $ is defined as equation (12). As usual, r c 
is the crossover scale in the DGP model and a measures the strength of the GB curvature 
effect on the brane. If a = 0, the model reduces to the DGP model with a non-minimally 
coupled scalar field on the brane, while for r c = we recover the Gauss-Bonnet braneworld 
model. Here, we restrict our study to the case where the bulk black hole mass vanishes and 
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we neglect the bulk cosmological constant too ( that is, A 5 = 0). For a spatially flat FRW 
brane, the Friedmann equation is given by 



3 V 2 



{H 2 - $) 



(42) 



We restrict ourselves to the normal branch of this GBIG-inspired scenario and therefore 
there is no ghost instabilities in this case. The energy density on the brane p, is composed 
of a CDM component ( with energy density p m ) and a cosmological constant A 4 



P = Pm + A 4 . 

Considering the normal branch of the equation (43), we have 

snG, . x a 4 mr'H s aH , 



H 2 = 



3 ipm + Pcp) + 



Comparing this equation with the following Friedmann equation 

rr9 8itG , . 8tvG , eff . 



we find 



8ttG e// ,_A 4 [/?($] 
3 (PdeJ " 3 



3 

- X H. 



(43) 



(44) 



(45) 



(46) 



where G = G eff = (8nM, 



r f 



and the crossover scale takes the following form 



2M| 



(3{4>)r c . 



(47) 



Equation (46) shows the dynamical screening of the brane cosmological constant due to both 
non- minimal coupling and the Gauss-Bonnet effect. In which follows, we set for simplicity 
P(4>) = |(1 — t,4> 2 ) where £ is the conformal coupling. We note also that in our numerical 

calculations we adopt the ansatz <p{i) = 0o^ M an d a(t) = (t 2 + jz^j , where p and v are 
positive constants. We obtain from equation (46) 



8nG {peM) = A, _ + 8 



+ 3«^ 2 )- 



(48) 



Figure 11 shows the variation of p e ^ E versus ^ in a constant time slice. The range of £ are 
chosen from [67-71] constraint by the recent observations. As this figure shows, by increasing 
the values of the nonminimal coupling, p e ^ E decreases in a fixed redshift( fixed time slice). 
Figure 12 shows the variation of the effective dark energy density versus the redshift. As 
this figure shows, p EE increases by decreasing the redshift and this is exactly the phantom- 
like behavior we are interested in. Note that this phantom-like effects is realized without 
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introducing any phantom matter on the brane and only screening of the brane cosmological 
constant due to non-minimal coupling and the curvature effect causes such an intriguing 
aspect. 

Figure 13 shows the variation of the effective dark energy density versus the redshift and 
the non-minimal coupling. We note that while the introduction of a phantom field requires 
the violation of the null energy condition, here this energy condition is respected ( at least 
in some subspaces of the model parameter space) since we have not included any phantom 
matter on the brane. In other words, since the phantom-like dynamics realized in this setup 
is gravitational ( the quintessence field introduced on the brane plays the role of standard 
matter on the brane), the null energy condition could be respected in some subspaces of the 
model parameter space. 



rho 




Figure 13: Variation of the effective dark energy versus the redshift and the non- minimal coupling. 



6 Conclusion 

In this paper we have analyzed in some detail a modified AD GP- Gauss-Bonnet model which 
corresponds to a 5D braneworld model where the bulk is a 5D Minkowski space-time. This 
model contains a Gauss-Bonnet term in the bulk and an induced gravity term on the 
brane which is modified in the spirit of f(R) gravity. We incorporated also a canonical 
(quintessence) scalar field non-minimally coupled to induced gravity on the brane. Our 
analysis is performed for the normal (or non self-accelerating) ghost-free branch which we 
have assumed to be filled by cold deark matter (CDM) and a cosmological constant. We 
have shown especially that how the brane accelerates at late-time in this normal branch. 
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An interesting feature of this model is the role played by the extra dimension: it induces a 
mimicry of a phantom-like behavior without introducing any matter that violates the null 
energy condition on the brane. This phantom-like behavior happens without any super- 
acceleration of the brane. We should emphasize that the Gauss-Bonnet effect which induces 
the ultraviolet corrections on the brane, modifies the phantom-like behavior at earlier times. 
Therefore, a consistent DGP brane model which essentially is an infra-red modification of 
the General Relativity, should have ultraviolet modifications as well and this is achieved by 
incorporation of the Gauss-Bonnet term in the bulk action. We have incorporated also the 
possible modification of the induced gravity on the brane. This is equivalent to introduce a 
canonical scalar field on the brane with General Relativity as the background gravitational 
theory. By adopting suitable and observationally viable ansatz for f(R, (J>) and the scale fac- 
tor, we have shown that this model realizes the phantom-like behavior in some subspaces of 
the model parameter space without violation of the null energy condition. The phantom-like 
prescription in our proposed model has no break down in some appropriate subspaces of the 
model parameter space and it realizes a smooth crossing of the phantom-divide line. 
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